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CN ! Abstract 

2 . If muon neutrinos produced in cosmological sources oscillate, neutrino telescopes 

I can have a chance to detect r- neutrinos. In contrast to f^'s the Earth is completely 

Q^' transparent for z^t-'s thanks to the short life time of r-leptons that are produced in 

^ . Uj- N charged current interactions, r-lepton decays in flight producing another Uj- 

Q \ (regeneration chain). Thus, z^t-'s cross the Earth without being absorbed, though 

' loosing energy both in regeneration processes and in neutral current interactions. 

Neutrinos of all flavors can be detected in deep underwater /ice detectors by means 
^ \ of Cerenkov light emitted by charged leptons produced in u interactions. Muon and 

' r-leptons have different energy loss features, which provide opportunities to identify 

^ r-events among the multitude of muons. Some signatures of r-leptons that can be 

firmly established and are background free have been proposed in literature, such 
as 'double bang' events. In this paper we present results of Monte Carlo simulations 
of r-neutrino propagation through the Earth accounting for neutrino interactions, 
r energy losses and r decays. Parameterizations for hard part and corrections to 
the soft part of the photonuclear cross-section (which contributes a major part to 
r energy losses) are presented. Different methods of r-lepton identification in large 
underwater /ice neutrino telescopes are discussed. Finally, we present a calculation 
of Uj- double bang event rates in km^ scale detectors. 
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1 Introduction 



^ij, ^ oscillations should lead to the proportion (py^ : (p^^ : (py^ = 1:1:1 
for neutrinos produced in cosmological sources that reach the Earth, though 
the flavor ratio at production in typical sources is expected to be (p^^ : (p^^ : 
(py^ = 1:2:0. Identification of UHE/EHE r-events in deep underwater /ice 
Cerenkov neutrino telescopes (UNTs) [1-6] would confirm oscillations already 
discovered at lower energies [7-9]. Measuring the ratio between cosmological 
i/^ and Ut fluxes one could also exclude or confirm some more exotic scenario 
[10,11], such as u decay, in which (p^^ ^ (f)y^. 

At energies Ey <1 PeV a general approach to discriminate rare neutrino events 
from the huge amount of atmospheric muons present also at kilometer wa- 
ter/ice depths is to select events from the lower hemisphere. These can be 
produced by neutrinos, the only known particle that can pass through the 
Earth with negligible absorption below these energies. Nevertheless, u cross- 
sections increase with energy. For muon neutrinos, absorption is considerable 
above 1 PeV, depending on the v zenith angle, hence on the path-length trans- 
versed in the Earth. On the other hand, i/^-'s generate r-leptons via charged 
current (CC) interactions Ut N ^ t N in the Earth. Being a short-lived par- 
ticle, r decays in fiight producing another and (in ^ 35% of the cases) 
secondary muon and electron neutrinos which are generated in decay modes 
T ev^v^ (5e=17.84%) and r \xv^,v^ (5^=17.36%) [12]. Thus, neutrinos 
of all flavors undergo a regeneration process and the Earth is transparent for 
z^r's of any energy. Nevertheless, they loose energy through the Earth due to 
N interactions (CC and NC) where the hadronic showers take away part of 
the energy. A minor part of the energy is lost also due to r-lepton propagation 
before decay. Thus, calculations for all flavor neutrino fluxes at detector loca- 
tion that result from propagation of an initial i/t flux through the Earth must 
necessarily take into account neutrino interaction properties, r-lepton energy 
losses and r-lepton decay. Such calculations are reported, e.g., in [13-20], both 
for monoenergetic i^t beams and for a variety of models for cosmological neu- 
trino spectra. 

r-neutrinos can be detected in UNTs identifying r leptons. To our knowledge 
the possibility of discriminating r's from muons through their different energy 
loss properties has not yet been analyzed. In this paper we discuss energy 
loss properties and speciflc r-induced signatures, such as 'double bang' (DB) 
or 'lollipop' [21] events, which should be affected by neghgible muon back- 
ground. In [22] the results of a calculation of the number of DB events in an 
UNT are published, but authors considered only down-going i/^-'s and ignored 
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those passing through the Earth to the detector from the lower hemisphere. 
In [14,23,24] the measurement of the ratio of up-going shower-hke and track- 
hke events is proposed, instead of an event- by-event identification of r- or 
//-events. As a matter of fact, r-leptons, as well as muons, produce events of 
both kinds in a UNT. i/g's produce shower-hke events, as well. In case Un os- 
cillate into Ur, the fraction of shower-like events is larger than what expected 
for no oscillations. 

In this paper we describe the results, prehminarily given in [25,26], of a MC 
simulation of r-neutrino propagation through the Earth. Both monoenergetic 
neutrino beams and spectra predicted by Protheroe [27] and Mannheim et 
al. [28] (which were not discussed in [14,15]) have been considered. We have 
taken into account corrections to photonuclear interaction cross-sections which 
play the main role in r-lepton energy losses, compared to pair production and 
bremsstrahlung. We used results pubhshed in [29] where i) the soft part of 
photonuclear cross-section based on the generalized vector dominance model 
originally published in [30,31] is corrected and ii) the hard part is developed 
in the QCD perturbative framework. The soft part corrections to photonu- 
clear interaction, not accounted for in [13-20] , have lead to an increase of the 
total T-lepton energy losses by 20-30% in UHE/EHE range. Also we have ana- 
lyzed the possibility to distinguish UHE/EHE r-leptons from muons in a UNT 
by different energy losses along tracks through production of electromagnetic 
and hadronic showers generated in bremsstrahlung, e"'"e~-pair production and 
photonuclear interaction. 

In Sec. 2 we describe the the tools used for this calculation with particular at- 
tention to r-lepton energy losses for photonuclear interaction; in Sec. 3 results 
on propagation of UHE/EHE r-neutrinos through the Earth are reported. In 
Sec. 4 we analyze the possibility to distinguish r-lepton from muon in a UNT 
thanks to their different energy loss properties and we give results on a cal- 
culation of the rate of DB events coming from both hemispheres in km"^ scale 
UNTs; in Sec. 5 we conclude. In Appendix A the corrected formula for muon 
and r-lepton photonuclear cross-section is given. 



2 Tools for Uj- propagation in the Ecirth 

A MC calculation has been developed to account for the following processes: 

• NC neutrino interactions that cause neutrino energy losses. 

• CC interactions. For z/g's and z/^'s it has been assumed that resulting elec- 
trons and muons are absorbed, while for z/,-'s r-leptons are followed. 

• r-lepton propagation through the Earth accounting for its energy losses. 

• r-lepton decay resulting in another i/^ appearance and in f»35% of the cases 
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also in z/g or z/^ production, i/^'s generated in r decays have been reprocessed 
again through this chain of processes. Thus, the 'regeneration chain' v^- 
r — i> z/^ — i> ■ ■ ■ has been simulated. 



^^T0 — '■ — NC interaction, 

[l yp_ ^ looses energy 




NC interaction 
i>T looses energy 

iv^t interaction,] 
res jl ting or /iA 
disappears 



y<f CC interaction, 
r appears 



Fig. 1. Schematic view of processes induced by high energy z^r's when propagating 
through the Earth. 

We have assumed the Earth composition made by standard rock {A=22, 
Z=ll) of variable density with the Earth density profile published in [32]. 
Processes undergone by r-neutrinos during propagation through the Earth 
are shown schematically in Fig. 1. 



2.1 Neutrino interactions 



To simulate u N interactions we have used a generator developed by P. Li- 
pari and F. Sartogo which is based on cross-sections described in [33]. It has 
been adapted for our use up to the high energy region we are interested in 
{El, > 10^ GeV) adjusting the efficiency of the rejection technique method 
described in [34]. Both for NC and CC interactions we have accounted only 
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Fig. 2. CC (solid lines) and NC (dashed lines) deep inelastic scattering cross-sections 
for Vr, i^T, i^/x; and (cross-sections for electron (anti-)neutrinos practically coincide 
with ones for muon (anti-)neutrinos) for isoscalar medium with Z/A = 0.5. Lines: 
cross-sections with structure functions CTEQ3-DIS (used in this work); markers: 
cross-sections with CTEQ4-DIS (from [42]). 

for deep inelastic scattering since in the considered energy range other chan- 
nels contribute negligibly. For what concerns deep inelastic scattering, the 
generator uses the electroweak standard formula for the inclusive differential 
cross-section and it is based on the LUND packages for hadronization [34- 
36]. Structure functions CTEQ3-DIS [37] taken from PDFLIB [38] suited for 
the high energy regime have been used ^ . More recent versions of CTEQ 
PDFs [39-41] produce differences on cross-sections at the level of up to 25% 
at lO^GeV (see Fig. 2). 



^ It should be noticed, however, that the CTEQ3 parton function extrapolation at 
small Bjorken implemented in PDFLIB, is more similar to the CTEQ6 one 

rather than to the CTEQ4 and CTEQ5. 
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2.2 Tau-lepton energy losses 



The MUM package (version 1.5) has been used to simulate r-lepton propaga- 
tion through matter. Comparing to [43] (where the first version of the package 
originally developed for the muon propagation was described) the package has 
been extended to treat r-leptons, accounting for their short life time and large 
probability to decay in flight. Besides, the newest corrections to photonuclear 
cross-section have been introduced. Formulas for cross-sections for e'^e~-pair 
production, bremsstrahlung, knock-on electron production and stopping for- 
mula for ionization implemented in MUM 1.5 can be found in [43] (Appendix 
A) ^ where they are given according to [50-60] . 

Main improvements in MUM 1.5 concern photonuclear interaction of leptons 
with nucleons, in which virtual photons are exchanged. We are interested in 
the diffractive region of the kinematic variable space (g^ < 0, transferred 
energy u large, xbj ~ /2v small). The photonuclear interaction has been 
treated as the absorption of the virtual photon by the nucleon and, using 
the optical theorem, it has been connected with the Compton scattering of a 
virtual photon, 7* -|- A?" — > 7* -|- A^. The Compton scattering in the diffractive 
region has been described by the vacuum exchange which, in turn, is modeled 
in QCD by the exchange of two or more gluons in a color singlet state. This is 
possible because, in the laboratory system, the interaction region has a large 
longitudinal size, and the photon develops an internal structure due to its 
coupling with quark fields. In the diffractive region, 7*A^-scattering dominates 
the Compton amplitude, while the contribution to it due to the photon bare 
component is smaller. 

It has been shown by HERA experiments [61-66] that the picture of the usual 
soft diffraction, soft pomeron exchange, does not work well if the center-of- 
mass energy of 7* and the target nucleon A^ is very large (i.e. if and u 
are in the diffractive region). For a description of the data in the framework 
of the simplest Regge-pole model, one needs at least two pomcrons: soft and 
hard ones. A two-component picture of the photonuclear interaction in the 
diffractive region arises very naturally due to inherent QCD properties (asymp- 
totic freedom, confinement, color transparency) and is the common feature of 
some recent quantitative models [29,44,45]. Generally, the photonuclear cross- 
section integrated over Q"^ — —q^ is expressed through the electromagnetic 



^ When dealing with r-lepton propagation one must change muon mass for r-lepton 
mass in all the formulas except for expression for = 1.9m^/Z^/^ in formula for 
bremsstrahlung cross-section (see page 074015-15, Appendix A in [43]) where the 
muon mass was introduced just as a mass-dimension scale factor. 
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structure functions of nuclei aL,T by the formula 




where Q^m — — v), v = AE/E is the fraction of energy lost in the 

interaction by a lepton of energy E with mass mi and a —1/137 is the fine 
structure constant. 

In the first version of the MUM code [43] only the soft part of this integral 
was accounted for according to Bezrukov-Bugaev (BB) parameterization that 
was developed in the frame of the generalized vector dominance model [30,31]. 
In MUM 1.5 two main improvements have been done. Firstly, new terms have 
been included in the BB parameterization [29]. Since these produce negligible 
efi^ects in the case of muon, they were not introduced in the original formula 
given in [30,31]. Nevertheless, these terms are essential in case of r-lepton 
(due to its larger mass compared to muon: m^u^O.lOSTGeV; mT-=1.777GeV) 
and they increase total r-lepton energy losses by ~20-^30% in UHE/EHE 
range. Secondly, the hard component of the photonuclear cross-section has 
been introduced using results obtained in [29]. This part is described by the 
phenomenological formula based on the color dipole model [46,47]. The main 
element of the approach used in [29] is the total cross-section a{r±, s) for 
scattering of dipolcs (gg-pairs) of a given transverse size r± from a target 
proton at fixed center-of-mass energy squared s. The r_L-dependence of cr[r±, s) 
is qualitatively predicted by perturbative QCD: 

a{r±,s) ~ rle ^(ris)^«. (2) 



In this formula Xh and r^^ are parameters determined from a comparison of 
the total photonuclear cross-section (or structure functions) with HERA data 
at small xbj- Such an expression for a{r±, s) was suggested in [48]. Besides, for 
a region of large s (outside of the HERA region) the unitarization procedure 
is used in [29] , which leads, asymptotically, to a logarithmic s-dependence of 
the dipole cross-section, 

(7""^*(ri,s) - Aj/lns. (3) 



The hard part of the photonuclear cross-section is shown in Fig. 4 and Fig. 5 
in [29] for muons and r's, respectively. We have made a polynomial parame- 
terization for this part that is given in Appendix A along with the complete 
formula for photonuclear cross-section used in this work. 
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Fig. 3. r-lepton energy losses in standard rock. Upper plot: energy losses due to 
photonuclear interaction. Lower plot: energy losses due to all radiative processes, 
including bremsstrahlung and direct e~'~e~-pair production computed according to 
[50-56]. In each plot: dashed line is for non corrected soft part [30,31]; dash-dotted 
line for corrected soft part [29-31]; dotted line includes hard component of photonu- 
clear cross-section according to [49]; solid line: includes hard component according 
to [29] (as it is treated in this work). 



r-lepton energy losses in standard rock {A=22, Z=ll) for photonuclear in- 
teraction and for the sum of all radiative processes are presented in Fig. 3. 
The importance of both corrections described above compared to the soft 
BB formula can be clearly seen. Results published in [49] for r-lepton energy 
losses including another calculation of the hard component of the photonuclear 
cross-section are shown for comparison, as well. The main difference between 
[29] and [49] is due to the fact that corrections to the BB formula (soft part) 
[30,31] have not been introduced in [49]. The hard component of the photonu- 
clear cross-section differs in [29] and [49], as well, since theoretical approaches 
applied in these two works are completely different. 



8 



2.3 Tau-lepton decay 



To generate r-lepton decays we used the TAUOLA package [67] which was 
developed for SLC/LEP experiments where r-leptons are produced in e'^e~ 
colhsions: e^e^ ^ Z ^ t^t~ . TAUOLA generates decays of r^'s of a given 
energy, taking into account all the effects of r spin polarization. 22 decay 
modes are treated (sum to almost 100% of the total width) including modes 
that are responsible for i/g and appearance: 

T~ e~ Pe (t"^ e+ (4) 

T~ jl" Vj. (r+ /X+ Ur) (5) 

In our simulation we tracked only neutrinos (of all flavors) resulting from r 
decay, since ranges of charged leptons are negligibly small compared to the 
Earth dimensions. Decay lengths of K^s and tt's resulting from r decay are 
much longer than their interaction lengths in the considered energy range, 
hence we did not simulate their decay neglecting secondary neutrinos ^ . 



3 Results on propagation through the Earth 

In this section results on the simulation of propagation of Ur and Ur monoen- 
ergetic beams through the Earth are presented where all processes shown in 
Fig. 1 have been accounted for. 

3.1 Results on monoenergetic tau-neutrino beams 

The outgoing spectra of initially monoenergetic z/^-'s after propagation through 
the Earth for three nadir angles of incidence on the Earth are shown in Fig. 4 
for two neutrino energies: 10 TeV and 1 EeV. Moreover, the outgoing spectra 
of secondary i/^'s and fg's that are produced in decays of r-leptons generated 
in i/r CC interactions are presented. Results for i/j- and beams are identical. 
At Eij^ =10^GeV neutrino cross-sections are small and most of r-ncutrinos 
cross the Earth without interacting keeping their initial energy (see the peaks 
in the upper panels in Fig. 4). The fraction of such neutrinos is lower at 
the nadir direction {9 = 0°) since the amount of matter crossed is maximum. 
This fraction increases with nadir angle. At E^,^ — lO^GeV all r-neutrinos 
undergo at least one interaction and the peak at initial energy disappears 

^ By the same reason we did not follow K^s and tt's produced in u CC interactions. 
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Fig. 4. Monoenergetic Vr beams incident on the Earth at nadir angles 0° (left), 
30° (center), 60° (right). Results after propagation through the Earth. Upper row: 

=10 TeV, lower row: = 1 EeV. Upper spectra: 1/^; lower spectra: secondary 
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Fig. 5. Left column: Mean number of secondary f^'s (= t'e's) per Ur that are gen- 
erated when a z/,- crosses the Earth undergoing regeneration processes vs E^^ (solid 
lines); mean number of secondary i^^'s (= fc'e's) that survive and come out from 
the Earth (dashed lines). Right column: mean number of CC (solid lines) and 
NC (dashed lines) interactions that occurs in the chain Uj- r v-j- - ■ ■ when 
a r-neutrino crosses the Earth vs its energy. Results for Uj- and v-r are identical. 

since neutrinos loose energy both in NC interactions and in regeneration chain. 
The fraction of secondary t/g's and t'^'s produced by the incoming t/^ is 5% 
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at =10^ GeV and 22% at E^^ =10^ GeV for the nadir direction ^ 
numbers are in good agreement with results pubhshed in [15]. 



These 




0.2 0.4 0.6 0.8 1 

D / D^^j.^^ (mwe / mwe) 

Fig. 6. Mean energy degradation for Vr with initial energies (from bottom to top): 
10^ GeV, 10^ GeV, 10*^ GeV, lO"^ GeV, 10^ GeV, 10^ GeV when penetrating the Earth 
at nadir angle 0° vs fraction of Earth diameter crossed. 

The fraction of secondary neutrinos that accompany r-neutrinos emerging 
from the Earth increases with initial r-neutrino energy but saturates at the 
level of 0.22 at some critical energy which depends upon nadir angle (Fig. 5, left 
column). These results are also in a good agreement with [15]. The CC cross- 
section increases with energy and, consequently, the number of secondaries 
generated along r-neutrino path increases also. Nevertheless, this growth is 
moderated by absorption of secondary i/^'s and z/g's which, in contrast to i/,-'s, 
do not regenerate. In the right column in Fig. 5 the mean number of NC and 
CC interactions occurring to a primary z/^ {vr) when it travels through the 
Earth is shown The number of CC interactions corresponds to the number of 
regeneration steps. This number grows both with r-neutrino energy and with 
amount of matter crossed, which is maximum at the nadir. 

Fig. 6 shows r-neutrino energy losses due to regeneration process and NC 
interactions when crossing the Earth. 



^ We consider all secondaries collinear with respect to the primary u^, an approxi- 
mation that at the energies to which we are interested is reasonable. 
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3.2 Results for some astrophysical diffuse v spectra 



We performed MC simulations for two spectra of cosmological neutrinos: the 
spectrum predicted by Protheroe [27] for an optically thick proton blazar 
model and an upper bound (not a source model) on diffuse v spectrum for 
optically thick sources [28] (MPR bound) ^ (see Fig. 7). 
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Fig. 7. Protheroe diffuse neutrino spectrum [27], MPR Umit on diffuse neutrino 
spectrum [28] and 'Quasar background flux prediction' (SS QC) [68]. 

We have considered neutrinos of all flavors in the incoming flux including z/g's 
and z/^'s produced in sources and i/,-'s that appear on the way to the Earth 
due to z/^ ^ oscillations. Neutrinos and anti-neutrinos of all flavors have 
been assumed to be present in the proportion 0^, : 0p = 1 : 1 in astrophysical 
neutrino fluxes. The background of z/^^e atmospheric neutrinos from pion and 
kaon decays (conventional neutrinos), as well as the contribution of prompt z/'s 
from charmed mesons (calculated in the Recombination Quark-Parton Model 



^ The 'Quasar background flux prediction' [68] was not considered, though it is 
more optimistic, since it is at the exclusion level by AMANDA [69,70] and Baikal 
[71] experiments. 
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frame) was simulated using spectra taken from [72]. 
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logi„[E„ (GeV)] log,„[E„ (GeV)] 



Fig. 8. Spectra of cosmological i^'s sampled according to Protheroe spectrum ([27]) 
incident on the Earth with 4 different nadir angles and after propagation through it. 
Left column: no oscillations, (p^e '■ '■ (pur = 1:2:0, right column: oscillations with 
maximal mixing, (j)^^ : (pi^^ : = 1:1:1. Thin solid lines: total incoming flux of as- 
trophysical i^'s; thick solid (right panels), dashed and dash-dotted lines (overlapping 
in right panels): out-coming i/^ -|- i/^, i'fj_ + u^, Ue + i^e after propagation through the 
Earth, respectively. z/^_e spectra include secondary neutrinos produced by i/^'s. In 
the left part of each panel atmospheric neutrino spectra (conventional -l-prompt [72] ) 
after propagation through the Earth are shown, upper thick markers: y^^ + v^, lower 
thin markers: + Ve- Atmospheric neutrino spectra are given up to ~300 TeV 
where cosmological neutrino fluxes start to dominate since the simulation statistics 
at higher energies is very poor due to the steepness of spectra. 

Results for Protheroe and MPR spectra for astrophysical neutrinos incident 
on the Earth with 4 different nadir angles and transformed after propagation 
through the Earth are shown in Fig. 8 and Fig. 9. r- neutrino fluxes exceed z/^^e 
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Fig. 9. The same as in Fig. 8 but astrophysical neutrinos are generated according 
to MPR upper limit on diffuse neutrino flux [28]. 

ones remarkably up to nadir angles 9 ~ 60° since i/^-'s are not absorbed by the 
Earth in contrast to z/^,e's. But their spectra are shifted to lower energies with 
respect to initial spectra due to energy degradation in regeneration processes. 
For all 6', the outgoing flux of astrophysical neutrinos exceeds the background 
of atmospheric neutrinos at E^^ > 10^ GeV. This cross-over determines the 
energy threshold for detection of diffuse neutrino fluxes in UNTs. Results for 
Protheroe spectrum are in a qualitative agreement with ones published in 
[13] (see Fig. 8 and Fig. 9 there) where different models for neutrino interac- 
tions and r-lepton energy losses were used. The fractions of secondary t'^^e's 
produced in regeneration chains with respect to the total i/^^g flux emerging 
after propagation through the Earth (made of of primary z/^,e's + secondaries) 
are 0.57, 0.18, 0.06, 0.02 (spectrum [27]) and 0.18, 0.06, 0.02, 0.01 (spectrum 
[28]) for nadir angles 6 =6°, 60°, 80°, and 85°, respectively. These numbers 
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Fig. 10. The ratio of Ur + (upper curves), + (middle curves) fluxes after 
propagation through the Earth and of secondary + i/^ (= + Pg) (lower curves) 
over incident flux oc E~'^ for four nadir angles. Markers: this work, lines: results 
published in [16] for nadir angles 0=0°, 9=30° and 6* =60° (data for 61=80° are 
absent in [16]). 

are larger than corresponding fractions for power-law spectra (pu oc and 
(pi, (X E~'^ reported in [16] ^ . On the other hand, the comparison of the results 
on attenuated neutrino fluxes obtained by our algorithm with the ones ob- 
tained in [16] for power-law spectra shows a reasonable agreement (Fig. 10). 
Hence, we can conclude that larger fractions of secondaries result from the 
fact that MPR and Protheroe spectra are harder compared to ones considered 
in [16]. 



" Authors do not provide numbers, nevertheless one can estimate them from Fig. 
1 in [16]. 
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4 Tau leptons: detection signatures 

4-1 Tau lepton identification through energy loss properties 

At energies larger than ^2-10^ GeV r-lepton ranges become comparable 
to the typical linear dimensions of operating and proposed UNTs (see Fig. 11) 
and r tracks can be reconstructed. As a first guess, one could try to distinguish 
a r-lepton track from a muon one by means of differences in energy losses. The 
larger r mass affects bremsstrahlung, e"'"e~-pair production and photonuclear 
processes differently respect to muons. As a result, e"'"e~-pair production dom- 
inates muon energy losses up to ^10^ GeV while r-lepton photonuclear 
interaction plays the main role above E^ ^ 10^ GeV (see Fig. 12) compared to 
other radiative processes. Hence, since probabilities of shower production for 
muons and r- leptons at a given energy are different, one could hope to select 
r-lepton events studying the distribution of showers produced in photonuclear 
or electromagnetic interactions along tracks even if the detector resolution is 
not so good to discriminate hadronic and electromagnetic showers. 



Fig. 11. Mean ranges for muons and r-leptons in standard rock and sea water. 
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Fig. 12. Muon (upper plot) and r-lepton (lower plot) energy losses in standard rock 
due to bremsstrahlung (dash-dotted lines), direct e"'"e~-pair production (dashed 
lines), photonuclear interaction (dotted lines). Thick solid lines stand for total ra- 
diative energy losses. Bremsstrahlung for r-lepton is strongly suppressed due to 
factor, therefore the bremsstrahlung curve almost coincides with the x-axis. 

Nevertheless, the results of the calculation that we have performed for r- 
leptons with energies in the range 10^ GeV< Er <10^ GeV indicate that such a 
method works badly because the largest differences between shapes of r-lepton 
and muon cross-sections lies in a range of relatively large v = AE/E, where 
interaction lengths exceed the detector size. The numerical data on r-leptons 
of energies 10^ GeV and 10^ GeV and muons of lower energies (1.10-10^ and 
1.56-10®, respectively) are presented in Fig. 13 and Fig. 14, which show values 
of energy losses and radiation lengths for different ranges of energy transfered 
to secondaries {lO^GeV < AE < lO^GeV, lO^GeV < AE < W^GeV, 
lO^GeV < AE < lO'^GeV). Muon energies were chosen so that total muon 
energy losses are equal to total energy losses of considered r-leptons. Data 
on electromagnetic processes (bremsstrahlung + e"'"e~-pair production), pho- 
tonuclear interaction and all interactions are presented separately. One can 
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Fig. 13. Left column: energy losses in sea water due to electromagnetic in- 
teractions (upper panel), photonuclear interactions (middle panel) and all in- 
teractions (lower panel) for energy transferred to secondaries belonging to dif- 
ferent ranges {lO'^GeV < AE < W^GeV, W^GeV < AE < W'^GeV, 
lO'^GeV < AE < WGeV). Filled boxes correspond to a r-lepton with E^ = 
W GeV, empty ones to a muon with energy E^ = 1.1 10^ GeV. Right column: 
Mean interaction length of muon/r-lepton in sea water for different ranges of en- 
ergy transfered to secondaries. The horizontal lines indicate 100 m and 1km which 
are the typical linear sizes of existing and planned UNTs. 



see that 'partial' energy losses and radiation lengths of r-leptons and muons 
of lower energies differ remarkably only for large transferred energies when 
radiation length exceeds 100 - 1000 m, typical linear dimensions of existing 
and planned detectors. Thus, there are too few high energy showers inside the 
detector sensitive volume to provide enough statistics to distinguish r-lepton 
and muon tracks. Even in the hypothesis that the detector capability is so 
good to distinguish electromagnetic showers from hadronic ones, it is very 
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Fig. 14. Tlie same as in Fig. 13 for a 10^ GeV r-lepton and a 1.56 • 10^ GeV muon. 

difficult and, most probably, impossible to distinguish a r-lepton track of a 
certain energy with respect to a muon track of ~6^11 times lower energy. This 
difficulty also concerns 1 km^-scale UNTs. Of course, the observation of a very 
high energy shower with reconstructed energy so large to be inconsistent with 
the energy reconstructed using the rest of the track could be an indication in 
favor of a r-lepton. Nevertheless, the probability of such a shower occurrence 
is too small that this cannot be considered a standard reliable method. 



4-2 Topological signatures 



As a matter of fact, signatures such as 'double bang' (DB) or 'lollipop' events 
proposed in [21] seem to be smoking guns to recognize r-lepton events in 
UNTs. DB event features are an hadronic cascade at the u N CC interaction 
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vertex and an hadronic or electromagnetic cascade corresponding to r-lepton 
decay. The r-lepton track connects both cascades ^ . A lollipop event con- 
sists of a partially contained r-lepton track with the first cascade at the v 
interaction vertex outside of the UNT sensitive volume and with the second 
cascade produced by r-lepton decay contained in the volume. On the other 
hand, the inverse pattern with only the first cascade visible can not be dis- 
tinguished from Vp,N interaction followed by the muon track. Both 
DB and lollipop signatures provide a clear and background-free evidence of 
Vj- detection. As a matter of fact, for r decaying into electronic and hadronic 
channels there is no muon track following the second cascade. Besides, DB 
events do not contain an in-coming track before the first cascade. In principle, 
an up-going muon generated in a CC interaction could mimic such a sig- 
nature if it would produce a shower that takes all the muon energy inside the 
instrumented volume. But the probability to have such an event for a muon 
track segment of the order of 1 km is less than 5-10""^ in water for energies 
larger than 1 PeV. This value corresponds to the probability that a muon 
looses more than 99% of its energy in one interaction. This number can be 
considered as an upper limit on the signal to noise ratio, where the signal 
are DB/loUipop events and the noise are up-going i/-induced muons. In fact, 
astrophysical fluxes are lower or at most equal to expected astrophysical Vj- 
ones after propagation in the Earth and atmospheric neutrino fiuxes are even 
lower above 1 PeV (see Fig. 8 and Fig. 9). Another background source could 
be due to atmospheric muons that are dominated by prompt muons in the 
considered energy range. DB events cannot be reproduced since in any case 
the muon track before the first cascade would be visible. Nevertheless, a muon 
can mimic a lollipop event if it enters the detector instrumented volume and 
undergoes an interaction with such a large energy transfer that the muon stops 
after shower production. We estimate a rate of ~ 50 km~^ yr"-*^ considering 
the probability upper limit of 5 • 10~"^ (corresponding to more than 99% muon 
energy loss) and the prompt muon fluxes published in [73,74]. However, this 
is an upper limit since it is not possible to compute exactly the probability 
that a muon looses 100% of its energy since cross section parameterizations 
close to f = lS.EjE = 1 do not describe perfectly well the real cross sections. 
Of course more quantitative conclusions on possible backgrounds are needed 
that should be based on detailed simulations including detector response and 



One should note that at energies where the r-lepton track connecting the 2 cas- 
cades can be reconstructed in a UNT (E't- > 1 PcV), the r is not a minimum ionizing 
track, as instead it is assumed in [21]. In the example considered in [21] (a 1.8 PeV 
r-lepton) about 5.5-10® Cerenkov photons in the wavelength range 350 nm < A < 
600 nm are emitted by a 90 m long track. This number is ~ 300 times larger than for 
a minimum ionizing track. Nonetheless it is still several orders of magnitude lower 
than the number of photons from both cascades (~4-10"'^^ and ~10^^ from first and 
second cascades, respectively), thus the signature of DB remains unchanged with 
this correction. 
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reconstruction algorithms. At least qualitatively, DB and lollipop events are 
potentially "background free" , as it has been stated in the paper that proposed 
the existence of these topologies [21]. 

In 17% of the cases, that is for the r decay channel into a muon (see eq. 5), the 
second cascade is absent. In this case the r-lepton track is followed by a muon 
track without any shower. This topology seems to be not-recognizable as a r- 
lepton signature in a UNT, although the amount of light produced by showers 
along the r track (hereafter 'brightness') differs essentially from that along a 
muon track of the same energy. As an example, let us consider a r-lepton of en- 
ergy Et-= 10^ GeV that decays into a muon. For a simple estimation, we assume 
that the muon takes 1/3 of the r energy, thus ^3.3-10^ GeV. From results 
in Sec. 4. 1 (see Fig. 13) we know that the brightness of a 10^ GeV r track is sim- 
ilar to that of a // track with = 1. 1 • 10^ GeV. After r-lepton decay, the track 
brightness increases by about a factor of ^3.3-10^ GeV/1.1 •10*' GeV=3. This 
factor depends on r-lepton energy and in the range 10^ GeV < < GeV 
it varies between 2 and 4. To recognize such a signature detector energy re- 
construction should be remarkably better than this factor, which is not always 
the case for neutrino telescopes (see for instance Fig. 4 from [75]). Thus, when 
computing DB or/and lollipop event rates in a UNT, results must be reduced 
by a factor -B_^= 0.83 to exclude the non-muonic branching ratio in r-lepton 
decay. 

Estimates of DB rates in km^ UNTs have already been presented in [22] for 
down-going z/^'s, while events from the lower hemisphere were not considered. 
Galculations of up-going z/^ propagation were done in [14]. Authors give a 
qualitative conclusion that though UHE/EHE i/^-'s are not absorbed by the 
Earth, their energies decrease; hence the expected amount of r events in the 
DB energy range is low. Calculations in [24] performed for a neutrino flux 
El dNy^ I dEy^ = 10~^ GeVcm"^ s~^ assuming z/ oscillations but without ac- 
counting for neutrino absorption in the Earth result in a rate of ^0.5 yr~^ both 
for DB and lollipop events in an IceCube-like km^ detector [4]. The measure- 
ment of the ratio between muon and shower events that may be an indirect 
signature of v-,- appearance was proposed in [14,24] as an alternative indirect 
way to detect r-neutrinos. Here we present results of a calculation of the DB 
rate for both the upper and lower hemispheres in a km^ scale UNTs using 
results reported in Sec. 3 for the spectra in [27,28] that were not considered 
in [14,22,24]. The rate of totally contained DB events in a UNT is given by: 

0(1) oo 

nDB = 2T:pNAB_^j J Veff{E,^,e) I{E,^,e) a^^'iE,^) dE,^d{cose) (6) 

Here Na is the Avogadro number, p is the medium density (we use p = 1 g cm~^ 
which is close to sea water/ice density), I{E^^, 9) is the differential v^. flux. The 
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Earth shadowing effect is accounted for in the case of the lower hemisphere (for 
which the upper hmit of the integral is 0, while cos = 1 is used to compute 
the number of events for the whole sphere). The effective volume is: 

Veff^S,{9){L-Rr{E,^)), (7) 



where Sp is the projected area for tracks generated isotropically in azimuth 
for a fixed 6 direction of the incident neutrino on a parallelepiped. The paral- 
lelepiped has the dimensions of an IceCube-like detector [4] (1x1x1 km^) or a 
NEMO-like one [5] (1.4x 1.4x0. 6km^). In eq. 7, L is the geometrical distance 
between entry and exit point from the parallelepiped, is the r-lepton range, 
(t'^^ is total CC u cross-section. Emin — '^ ■ lO^GeV corresponds to a r-lepton 
range of i?;"^" = 70m. 

Table 1 



Number of DB events in km^ detector per year. 



Spectrum 


IceCube-like (A/"_27r / / N4t,) 


NEMO-like {N^2n / / N^n) 


Protheroe 


0.6 / 1.3 / 1.9 


0.8 / 1.7 / 2.5 


MPR 


0.8 / 1.9 / 2.7 


1.2 / 2.6 / 3.8 



Table 1 shows DB rates for the lower, upper and for both hemispheres. Values 
for upper hemisphere are 3-r6 times lower compared to [22] since more opti- 
mistic predictions for diffuse neutrino fluxes [68] are used there. At the same 
time, our figures for DB rates in IceCube-hke UNT are 4-=-5 higher compared 
to results reported in [24]. There are at least four factors that contribute to 
this difference: 

• different u flux models are used; 

• minimum r range is i?^*" = 70 m in our calculations while in [24] it is equal 
to i?^*" = 200 -7- 400 m. Such a large minimum range was chosen according 
to the qualitative requirement to reconstruct the r-lepton track connecting 
the two cascades. This difference accounts for a factor of ~4 in the case of 
an E~'^ V flux between 10^ and 10^ GeV. Since horizontal photomultiplier 
spacing in IceCube detector is ~125m, but the vertical spacing is 17 m [4], 
we believe that this requirement could be too restrictive. The proper choice 
of the minimum range should indeed be done quantitatively using a full 
simulation of the detector including reconstruction procedures; 

• for the lower hemisphere we accounted for r-neutrino energy losses in the 
Earth. This losses can be such that events do not fall anymore above the DB 
'energy threshold' (that figures for lower hemisphere in Table 1 are about 
a factor of ^2 lower compared to the upper hemisphere ones). These losses 
were not accounted for in [24]; 

• in contrast to [24] we reduce all figures for rates by a factor B^^— 0.83 to 
exclude the muonic mode of r-lepton decay. 
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5 Conclusions 



Considering the AGN diffuse neutrino flux models in [27] and the upper limit 
on optically thick v sources in [28] we have found that in a km^ UNT one can 
expect 2-^4 DB events yr~^ ^ . This result depends on the assumed neutrino 
spectrum and on the minimum r-lepton range that can be reconstructed in the 
detector, hence on its geometry and reconstruction algorithms. Hence, these 
values should be considered as indicative, while full simulations for specific 
neutrino telescopes should be performed. Moreover, we have investigated the 
possibility of identification of r-lepton from muons by studying the distribu- 
tion of showers produced by electromagnetic and photonuclear interactions 
along the track. This seems to be very difficult. Detection of DB or lollipop 
events are the cleanest signatures even though rates seem to be quite small. An 
alternative indirect way to identify r's consists in measuring the ratio between 
shower and track events. In any case, simulations of r-neutrino propagation 
through the Earth should be performed including regeneration processes and 
r-lepton energy losses. For this, in particular, one needs to account for correc- 
tions to the soft part of the photonuclear interaction which increases the total 
r-lepton energy losses by Ri20-=-30% in UHE/EHE range. 
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A Appendix: photonuclear cross-section 

The differential cross-section for photonuclear interaction used in this work has 
been published in [29]. It is rewritten below in the form adopted in [43] (Ap- 
pendix A). In the formula the following values are used: a = 7.297353x 10^'^ - 
fine structure constant; A = atomic weight; mi = = 0.1057 GeV for muon 
and mi — m-r — 1.777 GeV for r-lepton; mjv = nucleon mass; tt = 3.141593. 



^ The rate of lollipop events is expected to be slightly larger than that of DB events. 
In fact, for lollipops the second cascade is not required to be inside an UNT, while 
for DB both should be contained. Hence, the detection probability for lollipop events 
is essentially larger than for DB at >10PeV [24], but due to the steepness of 
neutrino spectra the difference in rates is not very large. Nevertheless, an estimate 
would require a more detailed detector simulation than what has been done in this 
work. 
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Table A.l 

Coefficients for -^^^ = \ Y1\=q o,k logio formula for photonuclear 



cross-section. The upper figures stand for muons, lower ones does for r-leptons. 
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^(^) = 1--+^, 

V 



G{z) 



1 {l + z)e-^-l 

2 ^ i2 



G'(z)=3 (Z = l), 



= 0.00282A^/V. 



t = 



= 0.54 GeV^ = 1.80 GeV^ 



with total cross-section for absorption of a real photon of energy u = s/2mN ~ 
vE by a nucleon, a^N, parameterized according to [30,31] 

a^N = [114.3 + 1.6471n2(0.0213i/)]//6 



and hard part of cross-section described by polynomial parameterization 



(coefficients arc given in Tabic A.l). Parameterization for -^^^ is valid in 

the range lO^^ < < 1, 10^ < ^ < 10^ GeV , in MUM 1.5 ^^j^ is assigned to 
for £■ = 10^ GeV and v — 10~^, values between v — 10~^ and v — 10~^, as 
well as between E — 10^ GeV and E — 10^ GeV are computed by interpolation. 
Values for energies between ones that are given in the Table A.l are computed 
by interpolation, as well. 
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